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Abstract 

In this article, we take into account the one-pion exchange force besides the one- 
gluon exchange force to study the mass difference of the tt and p mesons with the 
Bethe-Salpeter equation. After projecting the Bethe-Salpeter equation into an simple 
form, we can see explicitly that the bound energy 1^ ^> \E P \. 
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The constituent quark models have given many successful descriptions of the hadron spec- 
troscopy, the simple constituent quark mass plus hyperfine spin-spin interaction model 
works well for the ground state mesons p], 

> 

o: Mm = Ml + M2 + C7 SS' (1) 

i> ■ 

where the coefficient C can be fitted phenomenologically. The masses of the ground state 
pseudoscalar and vector mesons are M n = 140 Me V, Mk = 494 Me V, = 548 Me V, 
M v > = 958 MeV, M p = 775 MeV, M w = 783 MeV, M K * = 892 MeV, and M = 1019 MeV 
[2]. The mass difference between the tt and p mesons are huge, we have to resort to large 
hyperfine spin-spin interactions for explanation. The fine spin-orbit interactions and the 
hyperfine spin-spin interactions are usually studied in the relativized quark model based 
on the one-gluon exchange plus linear confinement potential motivated by QCD [3], for 
more literatures, one can consult the comprehensive review Ref.[4]. One may wander 
why the contributions from the hyperfine interactions in Eq.(l) are so large, and how to 
understand them in the quantum field theory. 

At the energy scale p = 4:irf n « 1 GeV, the approximate chiral SUl(3) x SUr(3) 
symmetry is spontaneously broken to the SUv(3) symmetry by the small current quark 
masses m u , md and m a , and there appear eight Nambu-Goldstone bosons (in the following 
we will neglect the word Nambu for simplicity). The masses of the Goldstone bosons are 
related with the current quark masses through the Gell-Mann-Oakes-Renner relation [5]. 
The quark fields q(x) are usually decomposed as 

q(x) = exp [-i 75 ^(x)A a ] q{x) , (2) 

where the q{x) and £, a {x) denote the constituent quark fields (or Goldstone free fields) and 
the octet Goldstone boson fields, respectively [6], there exist interactions among the quarks 
and the Goldstone bosons. For example, the spectra of the nucleons, A resonances and 
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the strange hyperons are well described by the constituent quark model with the harmonic 
confinement potential plus the one-Goldstone-boson exchanges induced potential [7]. In 
the energy region between the confinement and the spontaneous chiral symmetry breaking, 
the elementary degrees of freedom are quarks, gluons and Goldstone bosons [8]. 

On the other hand, the mass breaking effects in the chiral doublets are very large [11 [9], 
for example, 

7r((T,140MeV) O / (600)(0+, 400 - 1200 MeV) , 
p(l",775MeV) o a x (1260) (1 + , 1230 MeV) , 

1+ 1- 
p(- ,938 MeV) <-> iV(1535)(- , 1525 - 1545 MeV) , (3) 

which requires that the chiral symmetry should be badly broken, here we use to denote 
the chiral rotations. The chiral massive quark dresses itself with the gluon cloud and 
quark- antiquark pairs, and acquires a dynamically generated large mass. We usually carry 
out the re-summation of the loops nonperturbatively with the Dyson-Schwinger equation, 
and obtain the Euclidean constituent quark masses by the definition p 2 = M 2 (p 2 ), which 
are compatible with the values used in the constituent quark models |10] [XT] . The light 
pseudoscalar mesons play a double role, as both Goldstone bosons and qq bound states. 

The exchanges of the one-gluon and one-Goldstone-boson between the two constituent 
quarks result in the hyperfine interactions He and Hp, respectively [7] fl~2]. 

Ha ~ „ , , r ^9 • A?<7j • a,- , 
MiMj 1 

and they both contribute to the spin-spin interactions. In this article, we take into account 
the contributions from the one-Goldstone-boson exchange force besides the one-gluon ex- 
change force, study the tt and p mass difference with the Bethe-Salpeter equation, and try 
to understand the difference in the quantum field theory. 

The Bethe-Salpeter equation is a conventional approach in dealing with the two-body 
relativistic bound state problems [13], and has given many successful descriptions of the 
hadron properties in a Poincare covariant way [111 114] . 

The article is arranged as follows: we solve the Bethe-Salpeter equation for the ud 
bound states in Sec. 2; in Sec. 3, we present the numerical results; and Sec. 4 is reserved for 
our conclusions. 
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2 Bethe-Salpeter equation 



We write down the ladder Bethe-Salpeter equation for the ud bound states in the Euclidean 
spacetimeU, 



S- 1 (q + CnP)x(q,P)Sf(q-CdP) 



d A k 



(27T) 



xa a, ^ Aa 9 2 s (q-k) 



2 ,MA - V ' ; 2 (g - k) 2 

aUq-k) 



' (q — k) 2 + m 2 



(5) 



Id 



MuH 
M u + M, 



the P p is the four-momentum of the center of mass of the ud bound states, the q p is the 
relative four-momentum between the u and d quarks, the x(g, P) is the Bethe-Salpeter 
amplitude of the ud bound states, and the g^iq — k) and g s (q — k) are the energy dependent 
7r-quark and gluon-quark coupling constants, respectively. In this article, we take the g 2 (k) 

as a modified Gaussian distribution, g 2 (k) = A (jz^j ex p( — Tv 7 )' wnere t ne strength 
A and the distribution width A 2 are free parameters. The ultraviolet behavior of the 
modified Gaussian distribution warrants that the integral in the Bethe-Salpeter equation 
is convergent. 

The Euclidean Bethe-Salpeter amplitudes of the ud bound states can be decomposed 

as 

X*(<1,P) = 7 5 {F w (q,P) + iPF?(q,P)+i<iq-PF£(q,P) + \P,d\F£(q,P)} , 
X p (q,P) = niF p (q,P)+PF 1 p (q,P)- ( iq-PFi(q,P)+i\P,d]F^q,P)} 
+q-e{q- PF p (q, P) + 2i PF p (q, P)} 

+q ■ e {F p (q, P) + i Pq ■ PF p (q, P) - i 4F p (q, P) + \P,d]F p (q, P)} , (6) 

due to Lorentz covariance [HUE], where the e p is the polarization vector of the p meson, 

the F^{q,P), F?{q,P), F p (q,P) and F p (q,P) are the components of the Bethe-Salpeter 

i 

amplitudes, which can be expanded in terms of Tchebychev polynomials Tn (cos 6) |16| . 

where 9 is the included angle between q p and P p . Numerical calculations indicate that 

i 

taking only the terms T 2 (cos#) = 1 can give satisfactory results [17]. If we take into 
account the small terms with n > 1, the predictions may be improved mildly. In the 
following, we use the amplitudes F n / p (q 2 , P 2 ) and F^ p (q 2 ,P 2 ) to denote the n = terms 

of the Bethe-Salpeter amplitudes F n / p (q, P) and F^ p (q, P), respectively. Then the Bethe- 
Salpeter equations can be projected into four and eight coupled integral equations for the tt 
and p mesons respectively, and it is very difficult to solve them numerically. Furthermore, 
we cannot obtain physical insight from those involved integral equations. 

Multiplying both sides of the Bethe-Salpeter equations of the tt and p mesons by 
75 y^P] [18], and f. \$, P] + \$, P] ft, q ■ eq ■ P P respectively, completing the trace in the 
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In this article, we use the metric 5 M „ = (1,1,1,1), {7 M 7„ + 7i,7 M } = 25 M „, the momentums fc M 



(&4, k), q M = (q 4 , it) and P M = (iE, ~P) with P 2 = —M, 



2 

Ip- 
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Dirac spinor space, carrying out the integrals for the included angle 9, and neglecting the 
small components Fg, F£, Fg and Fg, we can obtain the following three relations, 



FAQ 2 , P 2 ) - (M u + Mj) F?(q 2 , P 2 ) 
F p (q 2 , P 2 ) + (M u + M d ) FP(q 2 , P 2 ) 
2F p (q 2 , P 2 ) - (M u + Mj) F£(q 2 , P 2 ) 

The Bethe-Salpeter amplitudes can be approximated as 

iP 



0, 
0, 
0. 



(7) 



75 1 + 



M u + M d 
P 

M u + Mi 



fAq 2 ,p 2 ) 

2q ■ e 



+ 



M u + M d 



FJq 2 ,P 2 ), 



(8) 



if we also neglect the small components i^, F% and Fj . Then the involved Bethe-Salpeter 
equations can be projected into the following simple form, 



q 2 + M u M- d 



1 + 



F«{q 2 ,P 2 ) = 



d 4 k 
(2^ 



F n (k 2 ,P 2 



16 g 2 (q-k) , gl(q-k) 



+ 



q 2 + M u M- d 



(M u + M d ) 



F p (q 2 ,P 2 ) 



3 (q — k) 2 (q — k) 2 + m 2 
d '\F p (k 2 ,P 2 ) 



(2tt) 4 pv 

>g 2 (<i-k) gl(q-k) 



3 (q — k) 2 (q — k) 2 + m 2 



.(9) 



If we take q 2 = and = 0, and assume that there exists a physical solution, then 



M u M d 



M 2 , 

n/p 



(M u + M d ) 



/d 4 k 
—F n/p (k 2 ,-M 2 /p )G n/p (0 



k), 
(10) 



where the G n / p (k) denotes the interacting kernels. In numerical calculations, we observe 
that the Bethe-Salpeter amplitude F v / p (k 2 , —M 2 , ) has the same sign in the region k 2 > 0, 



Tr/p 



r d"i 
J W 



d 4 k F w/p (k 2 ,-M 2 ) 



(M u + M d ) 2 J (2Tr)*M u M d F n/p (0,-M 2 /p ) 
and obtain an simple relation (or constraint), 

M 2 /p <(M u + M d ) 2 , 



G n/p (0-k)>0, (11) 



(12) 



which survives for q 2 > (although the relation is not explicit for q 2 > 0), i.e. the 
bound energy E w / p originates from the interacting kernel G 7T / p (k) and should be negative, 
E n / p = M n / p — M u — Mj < 0. The numerical calculations indicate that above arguments 
survive in the case g 2 7^ 0. 
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From Eq.(9), we can see explicitly that the one-gluon exchange force in the tt channel 
is more attractive than that in the p channel due to the factors ^ and |; furthermore, 
the one-pion exchange force is attractive in the tt channel and repulsive in the p channel, 
the bound energy \E n \ S> \E p \ can be accounted for naturally. 

We can introduce a parameter A(P 2 ) and solve above equations as an eigenvalue prob- 
lem, the masses of the tt and p mesons can be determined by the condition X(P 2 = 



q 2 + M u Mg 



q 2 + M u Mi 



1 + 



1 + 



(M u + M- d 



F n (q 2 ,P 2 ) = X(P 2 



W) 4 



Fn(k 2 ,P 2 ) 



16 g 2 s {q-k) , gl{q-k) 



+ 



FJq 2 ,P 2 ) = 



3 (q — k) 2 (q — k) 2 + m 2 
9 2 {q-k) 9 2 M-k) 



3 (q-k) 2 (q-k) 2 +m 2 



(13) 



3 Numerical results 

The constituent quark masses of the u and d quarks are taken as M u = Mg = 400 MeV. 
The strength parameter A and the distribution width A are free parameters, we take the 
values A = 200 MeV and A = 146(105) for the one-gluon (one-pion) exchange. Other 
values of the A and A also work, we choose the present parameters for illustration. 

We solve the Bethe-Salpeter equations as an eigenvalue problem numerically by direct 
iterations, and observe the convergent behaviors are very good. The numerical results 
for the Bethe-Salpeter amplitudes are shown in Fig.l. From the figure, we can see that 
the Bethe-Salpeter amplitudes center around zero momentum and extend to the energy 
scale about q = 0.4 GeV and 0.7 GeV for the p and tt mesons respectively, the stronger 
interactions in the tt channel result in more stable bound state than that in the p channel, 
as the bound energies are E n = —660 MeV and E p = —25 MeV, respectively. In numerical 
calculations, we observe that the one-pion exchange force plays an important role and 
should be taken into account. 

With the following simple replacements in Eq.(13), 

M j ->• M- s = 536 MeV , 

M w -»• M E - = 494MeV, (14) 

where the A = A ^ M m + +m* ) denotes the SUv(3) breaking effect for the coupling con- 
stant, we can obtain the corresponding solutions for the pseudoscalar meson K and vec- 
tor meson K* with the eigenvalues \{P 2 = -M\ = —(494 MeV) 2 ) = 1 and \{P 2 = 
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Figure 1: The Bethe-Salpeter amplitudes of the bound states. 

-M 2 K , = -(892 MeV) 2 ) = 1, respectively. The bound energies are Ex = —442 MeV and 
Ek* = —44 MeV, respectively, which indicate that there exists a more stable bound state 
in the pseudoscalar channel than that in the vector channel. 

4 Conclusion 

In this article, we take into account the one-pion exchange force besides the one-gluon 
exchange force to study the mass difference of the tt and p mesons with the Bethe-Salpeter 
equation. After simplifying the involved Bethe-Salpeter equations, we observe that the 
one-gluon exchange force in the tt channel is more attractive than that in the p channel, 
while the one-pion exchange force is attractive in the tt channel and repulsive in the p 
channel, the bound energy \E n \ » \E p \ can be accounted for naturally. 
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